Critical fields for vortex expulsion from narrow superconducting strips 
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We calculate the critical magnetic fields for vortex expulsion for an infinitely long superconducting 
strip, using the Ginzburg-Landau formalism. Two critical fields can be denned associated with the 
disappearance of either the energetic stability or metastability of vortices in the center of the strip 
for decreasing magnetic fields. We compare the theoretical predictions for the critical fields in the 
London formalism with ours and with recently published experimental results. As expected, for 
narrow strips our results reproduce better the experimental findings. 
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I. INTRODUCTION 

The Meissner effect [3], i.e., the complete expulsion 
of an external magnetic field when a superconductor is 
cooled below the critical temperature, T c , is a fundamen- 
tal property of the superconductivity. The ratio between 
the superconducting characteristic lengths, A and £ (pen- 
etration and coherence length, respectively), allows for a 
classification of the superconductors as type-I or type-II. 
The type-I ones (k = A/£ < 1/V2) exhibit a complete 
Meissner effect below a critical field. The existence of im- 
purities, grain boundaries, or simply geometrical effects 
can, nevertheless, favor the penetration of the magnetic 
field inside a type-I sample. Instead, type-II supercon- 
ductors allow intrinsically for the penetration of magnetic 
field filaments with quantized flux (vortices) which form 
a triangular lattice [2j. 

The behavior of a vortex in a superconducting film 
was first described by Pearl Q in the London limit. He 
found the penetration length for a thin film with thick- 
ness d to be described by A = X 2 /d. Also using the 
London approach for an ideal semi-infinite superconduc- 
tor Bean and Livingston^ first showed that a surface 
barrier is produced as as result of the competition be- 
tween the attraction of the vortex towards the surface 
due to an image antivortex outside the superconductor 
and the repulsion of the vortex due to surface screening 
currents that circulate in opposite direction to the vortex 
currents. Geometrical barriers for thin flat superconduc- 
tors were introduced trying to explain the magnetization 
of high critical temperature superconductors [5|. In the 
London limit various other problems such as the behav- 
ior of the vortex lattice in a type-II superconductor strip 
with current @ , the energetics of a vortex near the edge of 
a thin superconducting film at zero field Q , and the study 
of the Bean-Livingston barrier[8[ have been considered. 

Whether or not vortices can be found inside a type- 
II superconductor and the density of them depends on 
whether the superconductor is field cooled or zero field 
cooled before applying the magnetic field. The fact that 
the number of vortices is not a unique number as a func- 
tion of temperature and field is related, particularly in 



mesoscopic superconductors, to the existence of surface 
barriers for vortex escape and entrance 0, [lo| . The use 
of superconducting strips is common in some technolog- 
ical applications. Therefore, the characterization of this 
type of geometry is necessary for the design of supercon- 
ducting devices. The noise generated by the movement of 
vortices causes loss of coherence in qubits and also affects 
the sensitivity of devices as the SQUIDS[12||. In a recent 
work, Stan et al.\l^ have measured the minimum field 
below which the vortices are completely expelled from 
bidimcnsional Nb superconducting strips. The strips are 
cooled down through the critical temperature in the pres- 
ence of a magnetic field applied previously. When the 
magnetic field is greater than B m , vortices can be ob- 
served by a Scanning Microscopy Hall probe. For lower 
magnetic fields, vortices are completely expelled. 

From a theoretical point of view it makes sense to de- 
fine two critical fields which are associated with the disap- 
pearance of both the energetic stability and metastability 
of vortices in the center of the strip for decreasing mag- 
netic fields. While an upper critical field can be obtained 
simply in terms of ground-state considerations, a careful 
computation of the surface energy barrier in the deter- 
mination of the lower critical field is mandatory. In this 
work we study the nature of these two fields in narrow 
strips with the help of Ginzburg-Landau theory. We first 
review previous work based on the London model in Sec. 
II. In Sec. Ill we review the basics of Ginzburg-Landau 
theory applied to superconducting strips. In Sec. IV we 
make use of the latter approach and compare the theo- 
retical predictions for the critical fields obtained in the 
London formalism with ours and with recently published 
experimental results. As expected, for narrow strips our 
results reproduce better the experimental findings. The 
conclusions are presented in Section V. 



II. CRITICAL FIELDS IN THE LONDON LIMIT 

Following the Bean-Livingston model[3], the free en- 
ergy for a vortex in a strip as a function of the transversal 
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position x can be described by ll|, [l 
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where /io is the permittivity of vacuum, $o is the flux 
quantum, B is the applied field, and w is the width of 
the strip. Similar expressions for the vortex energy in 
a strip have been derived in Refs. 0, S The first term 
corresponds to the nucleation energy of the vortex and 
the second to the interaction with the screening currents. 
Equation Q] diverges for x = and is not valid within a 
distance £ from the edges. This is not surprising since 
it is based on the London approximation which neglects 
the core of the vortex. From the previous equation two 
critical fields, Bq and B s , can be defined. For low mag- 
netic fields the free energy G(x) displays a maximum in 
the center of the strip. As the magnetic field increases 
the maximum disappears at a field Bq 
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where G(x) presents a plateau around x = 0. For B > Bo 
a local minimum appears in x — and a vortex can be 
placed there in a metastable situation. At a magnetic 
field B s > B Q 
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the minimum becomes an absolute minimum of energy. 
Above this field the free energy of the strip with a vortex 
is lower than the free energy without vortices. The fields 
Bo and B s are two possible critical fields for the expul- 
sion of the vortices. Arguments have been put forward 
in favor of Bp as the relevant field for the total vortex 
exclusion (ll, 14 1 and in favor of B a 

Experimental details on the manufacture and the set- 
up to observe the vortices in strips as well as measure- 
ments of B m are described in Ref. [l3|. The strips had a 
thickness of d =210 nm, a length of L y =4 mm, and a 
width of w =1.6, 10 and 100 /im, with a critical temper- 
ature, T c =8.848 K. In order to find the values of B m for 
every strip, Stan et al. count the number of vortices for 
different values of the magnetic field after cooling down 
the samples through T c . Their results show a linear de- 
pendence between the number of vortices and the mag- 
netic field for higher fields. For low values of the field 
the dependence of the number of vortices with the field 
deviates from linearity until the vortices are completely 
expelled from the strip. The value of B m is found by 
extrapolating the linear zone. The working temperature 
Tf is mantained close to T c where the mobility of the 
vortices allows them to explore the lowest possible en- 
ergy configurations without getting pinned by impurities. 
This temperature is T f = 8M5K. Taking £ =38.9 nm 
and k =5.0, £(T» =320 nm and H c2 (T f ) = 3.217 mT. 



-1.20 




0.35 0.36 0.37 0.38 

B (H J 



0.39 



FIG. 1: Ginzburg-Landau functional for a superconducting 
strip of width 5 £. The transition of the ground state from 
Meissner (solid line) to 1 vortex row state (dashed line) occurs 
at 0.372 H c2 . 



Thus, the narrowest strip has an approximate width of 
5f . For the strips with width 10 /im, the experimental 
results of B m and the theoretical prediction of B s based 
on Eq. Q]are basically in agreement with each other. For 
the widest strip A <C w, the theoretical B s must be cor- 
rected by a factor of three to account for the partial field 
screening. However, for the narrowest strip, the Lon- 
don theory badly fails. For example, the theory predicts 
Bq > B s when the strips are narrower than 5£. For nar- 
row strips with widths of the order of £ the finite size of 
the vortex core cannot be ignored which is likely to be 
the reason behind the failure of the theory and the dis- 
agreement with the experimental results. For this reason 
we propose to describe the narrow strips using Ginzburg- 
Landau theory. 



III. GINZBURG-LANDAU THEORY IN THE 
LOWEST LANDAU LEVEL APPROXIMATION 

The superconducting density is found by minimizing 
the Ginzburg-Landau functional [lij 
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with AG = G s — G n where G s and G„ are the Gibbs 
free energies for the superconductor and normal state, 

respectively, and n = (—ihV ~ ^-Aj is the kinetic 
operator. When the magnetic field is applied perpendic- 
ular to the superconducting thin film (situated on the 
xy plane) , the effective penetration length A can become 
greater than the width of the sample w, particularly close 
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to T c . In this case the screening of the magnetic field is 
weak and the field penetrates almost uniformly so we as- 
sume b = H = Bz where B is spatially uniform. From 

now on we use the magnetic length I = yj and the 

cyclotron energy E c = as length and energy units. 
The cyclotron frequency, is to c = e*B/m*c. 

We now expand the order parameter as a linear com- 
bination of eigcnfunctions, cf>k, of the kinetic operator 
with eigenenergies eu- Choosing the Landau's gauge for 
the vector potential A = Bxy, one finds the following 
expression for the eigenfunctions: 



4> k =Ae iky Xk(x), 
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where fc is the set of wavevectors in the direction y, Xki x ) 
are functions without nodes, and A is a normalization 
factor. These functions are calculated numerically and 
obey the usual boundary condition for a superconductor- 
vacuum system. In the lowest Landau level [TtJ the free- 
energy functional becomes 
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with AG expressed in units of 2tt-J^E c 



and E c = *%*, p 
in units of E C V. L y is the length of the strip (L y — > oo) 
and N c is the number of terms in the linear combination, 
at = a + efc is the condensation energy of the term fc in 
the expansion. Details of the minimizing process can be 
found in Ref. Ill 



IV. CRITICAL FIELDS IN THE 
GINZBURG-LANDAU FORMALISM 

We determine B s by finding the magnetic field at which 
the Gibbs free energies of a system with N c = 2 (one 
vortex row) and N c = 1 (Meissner state) are equal. After 
minimizing the Ginzburg-Landau functional, we found 
for the narrowest strip a value of B s =0.372i7 C 2 as shown 
in Fig. [I] For higher magnetic fields the penetration of 
another vortex row in the strip is not possible. For N c = 
2 (one vortex row) the diference between wavevectors in 
the expansion, Afc = fe2 — kl, establishes the physical 
distance between vortices in the row, Ay, 
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For example the wavevectors at B s are k = ztO.8151^ 1 
and therefore the distance between vortices is 3.85Z. The 
difference between wavevectors changes with the mag- 
netic field, therefore we have to search for the wavevectors 
that minimize the functional for every magnetic field. 
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FIG. 2: One vortex row energy as a function of the vortex po- 
sition across the strip where the energy of the Meissner state 
has been substracted. The magnetic field is Bo =0.265-ff C 2 
and the width 5£. The position is given in magnetic lengths. 
The insets show the superconducting density when the vor- 
tices are at the center and on the edge. 



By varying the values of the wavevectors, but keeping 
Afc fixed, we can calculate the free energy of one vor- 
tex row with respect to the Meissner state, G(x), as a 
function of the vortex row position across the strip. The 
variation of the wavevectors keeping their difference con- 
stant implies the numerical variation of coefficients and 
the minimization process turns to be completely numer- 
ical. In general for any configuration the coefficients of 
the linear combination are different between them and 
they are equal only when the vortex row is in the center 
of the strip. The displacement of a vortex row across the 
strip is associated with a current induced in the strip. In 
the Ginzburg-Landau formalism the existence of current 
in a superconducting strip with a vortex row increases 
the vortex linear density [191] and the distance between 
wavevectors in the expansion changes. Here we neglect 
this effect, which is small anyway. 

We determine the critical field Bo calculating G(x) for 
various magnetic fields. Bq is determined by the appear- 
ance of a plateau in the free energy in the center of the 
strip or the reversal of the sign of the second derivative 
d 2 G(x)/dx 2 \\ x=0 . We find B =0.265i7 c2 (see Fig. [2J. 
At Bo the distance between vortices is 4.269Z (see in- 
set). As expected, this distance is greater than that at 
B s due to the increasing vortex density at higher fields. 
At B = H c2 the distance between vortices is reduced to 
1.872/, but have not disappeared yet due to the presence 
of the surfaces. For magnetic fields above H C 2 the vor- 
tices merge creating a channel without superconductivity 
in the center of the strip 18]. We also compute G{x) for 
the critical field B s (see Fig. As expected a zero- 
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FIG. 3: One vortex row energy as a function of the vor- 
tex row position across the strip where the energy of the 
Meissner state has been substracted. The magnetic field is 
B s —0.372H C 2 and the width 5£. The position is given in 
magnetic lengths. The insets show the superconducting den- 
sity when the vortices are at the center and on the edge. 



energy mimimum is obtained in the center of the strip 
for which the superconducting density is shown in the 
inset. 

The free energies curves at both fields Bq and B s show 
a finite value when the vortices are placed on the edge. 
This behavior is different from that predicted by Eq. [T] 
and can be attributed to a correct description of the order 
parameter in the Ginzburg-Landau calculation. The free 
energy for a vortex on the edge is superior to the one 
of the Meissner state in both cases. The position of the 
vortex is associated with the position of the center of its 
core although the core can be partially inside the strip. 
For instance, when the center of the core is exactly placed 
on the edge of the strip, the vortex has half of the vortex 
inside the strip and the other half outside (see insets in 
Figs. [2] and [3}. We note finally that the free energy as a 
function of the position of the vortex presents a maximum 
near the surface for B s similar to the one obtained in the 
London approximation [il S [lO] . If a vortex nucleates at 
the surface still has to surpass the small barrier before it 
slides down towards the center of the strip. 

We now repeat the calculation of the critical fields Bo 
and B s for widths between 3£ and 5£. We find the ex- 
pected behavior of the critical fields for all the widths, 
namely B s > Bo, as shown in the Fig. dJ The closest 
result to the experimental data for the 5£ strip corre- 
sponds to the critical field B s using the formalism of 
Ginzburg-Landau. This result is compatible^ with the 
field cooled measurements of Stan et al. 



until the critical field Bq is reached due to the saddle 
point energy barrier between the vortex solution and the 
Meissner solution. This barrier avoids the transition be- 
tween the vortex states and the Meissner state despite of 
the former being energetically favorable. Thus, although 
the ground state for magnetic fields lower than B s has 
no vortices, the barrier would keep the vortices inside the 
strip until the critical field Bq is reached. A combination 
of other effects not considered in this work such as ther- 
mal fluctuations and defects on the edge would decrease 
the escape barrier and would favor B s over Bq also in 
decreasing fields. 

An almost linear behavior for the density of vortices as 
a function of the magnetic field can be seen in Fig. [5] for 
large values of B. We also find an almost linear behav- 
ior for the density of vortices as function of the magnetic 
field for wider strips and large fields with the exception of 
small steps produced by the entrance of individual vortex 
rows (not shown). The ground state (solid line) has no 
vortices for magnetic fields lower than B s although, as 
stated in previous paragraph, the system could evolve in 
a metastable state (dashed line) until Bo is reached. The 
arrows signal the region used by Stan et al. in their esti- 
mation of B m . For lower fields we observe a deviation 
from linearity which also seems to be the case in the Stan 
et al. experiment, although the authors attribute this to 
imperfections or pinning. We find a density of vortices 
greater than the experimental data, however the number 
of vortices expelled from the strip in the experiment is in 
agreement with the number the vortices expelled in our 
calculation in the same range of fields. The differences 
could be attributed to impurities or defects on the edges 
which are not accounted for in our simple geometry. 

Finally we would like to note that the results for the 
critical fields present a similar behavior to the one pre- 
dicted by London theory. According to London the- 
ory the expulsion field, B m , can be approximated to 
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where A; is a constant with a value close 



13fl. Never- 



theless, one should not forget that if the magnetic field 
decreases below B s the vortices would stay in the strip 



to <&o. In our case the same expression applies, but with 
an exponent for the denominator of 2.216. 



V. CONCLUSIONS AND DISCUSSION 

Two expulsion critical fields, Bq and B s , below which 
metastable and stable vortices, respectively, cease to exist 
have been determined for infinitely long superconducting 
strips and widths between 3£ and 5£ using Ginzburg- 
Landau theory. The expected behavior for the values 
of the fields is obtained, always being Bo < B s for this 
range of widths. The profile across the strip of the free 
energy of a vortex for the field Bo presents a flat zone 
in the center of the strip whereas the one for B s has a 
zero-energy minimum for vortices at the center. The the- 
oretical value of the two critical fields are closer to the 
experimental one than those obtained in the London for- 
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FIG. 4: Critical fields, B a and Bo, in the Ginzburg-Landau 
and London formalisms for different widths of a superconduct- 
ing strip. The theoretical result of B a in the Ginzburg-Landau 
formalism is the closest value to the experimental value (solid 
square) . 



by Stan et al. 13[ for wider strips. An overall agreement 
between our results for the number of vortices as a func- 
tion of the field and the experiment is obtained, including 
the low field deviation from linearity at small values. If 
the observed deviation is not due to imperfections or de- 
fects, the experimental value for the critical field would 
be smaller and in closer agreement with our calculations. 
Finally, one should keep in mind that, in perfect samples 
free of imperfections, and ignoring temperature effects, 
the critical field for the expulsion of the vortices should 
correspond to Bq, at least in decreasing fields. 
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FIG. 5: Linear density of vortices for a 5£ superconducting 
strip. Solid and dashed lines correspond to the stable and 
metastable one vortex row state, respectively. The system can 
sustain vortices inside the strip for magnetic fields lower than 
B s . Arrows delimit the field range used for the extrapolation 
of B m in Ref. [H 



malism. Our value for B s for a strip with a width of 5£ 
is the closest value to the experimental one, as suggested 
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